On the leptonic partial widths of the excited tjj states 
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The resonance parameters of the excited -(/j-faniily resonances, namely the ■!/;(4040), ipiAlQQ), and 
-0(4415), were determined by fitting the i?-values measured by experiments. It is found that the 
previously reported leptonic partial widths of these states were merely one possible solutions among 
a four-fold ambiguity. By fitting the most precise experimental data on the 7?-values measured by 
the BES collaboration, this work presents all four sets of solutions. These results may affect the 
interpretation of the charmonium and charmonium-like states above 4 GeV/c . 

PACS numbers: 12.38.Qk, 13. 66. Be, 14.40.Gx 
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I. INTRODUCTION 

Chamionium, bound state of a charm and an anti- 
charm quarks, is one of the most interesting two- 
body systems which were studied extensively in particle 
physics. Although the first charmonium state was dis- 
covered more than thirty-five years ago, there are still 
many puzzles in charmonium physics. The charmonium 
spectroscopy below the open charm threshold has been 
well measured and agrees with the theoretical expecta- 
tions (such as potential models and lattice QCD); how- 
ever, for the charmonium states above the open charm 
threshold, there are still lack of adequate experimental in- 
formation and solid theoretical inductions. For example, 
recently many new particles have been discovered, named 
XyZ-particles, and the overwhelming vector states in 
the 4 GeV/c^ to 5 GeV/c^ mass range make the classifi- 
cation of these states as the charmonia questionable 1- 
|3|. In explaining these vector charmonium states, the 
leptonic partial widths provide very important informa- 
tion. As we know, the vector quarkonium states could 
be either S-wave or Z?-wave spin-triplet states, with the 
5"- wave states couple strongly to lepton pair while the D- 
wave states couple weakly since the latter are only pro- 
portional to the second derivative of the wave-function 
at the origin squared, as expected in the potential mod- 
els. This leads people believe that the V'(4040) is the 35* 
charmonium state, ^(4160) the 2D state, and ?/'(4415) 
the 45* state. This has been a well accepted picture 
for more than two decades before the discovery of the 
so-called Y particles, namely, the y(4008) and F(4260) 
observed in e+e^ — > tt+tt^ J/f/; final state [J, Q, and 
the r(4360) and y(4660) observed in e+e" -^ n+TT-ip' 
final state [y, 0|- With seven states observed between 
4.0 GeV/c^ and 4.7 GeV/c^, some people started to cat- 
egorize some of these as non-conventional quarkonium 
states, while others tried to accommodate all of them in 
modified potential models. Many of the theoretical mod- 



els use the leptonic partial widths of these states to dis- 
tinguish them between S- and D-wave assumptions [8i,l9[, 
and most of the time, the values on the leptonic partial 
widths are cited from the PDG [1^ directly. Although 
the resonance parameters of these excited tp states have 
been measured by many experimental groups, all of them 
were obtained by fitting the i?- values measured in the rel- 
evant energy region. The most recent ones, which were 
from a sophisticated fit to the most precise i?- values mea- 
sured by the BES collaboration [ll|, ll2[, are the only 
source of the leptonic partial widths of these three tp 
states now quoted by the PDG [W^ . 

In fitting to the BES data, unlike the previous anal- 
yses, the BES collaboration considered the interference 
between the three resonances decaying into the same fi- 
nal modes, and introduced a free relative phase for the 
amplitude of each resonance |13j . The new parametriza- 
tion of the hadronic cross section results in a prononce 
increase of the -0(4160) mass, and significant decrease of 
the leptonic partial widths of 0(4160) and 0(4415). 

As has been pointed out in a recent study [14|, there 
are multiple solutions in fitting one dimensional distri- 
bution with the coherent sum of several amplitudes and 
free relative phase between them. Exactly the same kind 
of fit to the i?-value is performed in the present study, 
multiple solutions are indeed found in extracting the res- 
onance parameters of the excited states. The only dif- 
ference between these multiple solutions is the coupling 
to the e~^e~ , namely the leptonic partial width, while the 
masses and the widths of the resonances remain the same 
for all the solutions. 

In the following, we firstly introduce a simplified fit 
scheme similar to that used by the BES collaboration, 
and extract the resonance parameters. Then we study 
the multiple solution problem in the light of the toy sim- 
ulated data for illustrative purpose. We will discuss the 
consequence of the multiple solutions in fitting the R- 
values at the end of the paper. 
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II. THE FIT TO THE R DATA 

To faciliate our study, only the i?-values provided by 
the BES collaboration are used [O, [131 , and in the data 



fitting, only statistical uncertainties are considered, as 
the systematic errors at all the energy points are highly 
correlated. 

We fit the e'^e^ annihilation cross section cr(e+e^) = 
R ■ 86.85/s (s in GeV^ and a in nb). The standard 
chisquare estimator is constructed as 
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where a'^'^^' and Atr^^^^' indicate respectively the experi- 
mentally measured cross section and its error at the l-ih 
energy point (the number of points is denoted as iVpt), 
while a\^'^' is the corresponding theoretical expectation, 
which is composed of two parts 

a*''^-(s)=a'''='*-(s)+a™"-(s) . (2) 

Here cr'^""- is the contribution from continuum and is pa- 
rameterized simply as 



a™"-(s)=A + B(yj-2Afc±; 



(3) 



where A and B are free parameters, and M^ji is the 
mass of the charged D meson, cr'"'^*- is the contribu- 
tion from the resonances. Here following previous analy- 
ses |13l. Il5j. the assumption that the continuum produc- 
tion and the resonance decays don't interfere with each 
other is adopted. The three wide resonances shown in 
the data are close and have the same decay modes, the 
interferences between them must be included. Then the 
amplitude reads 



TA^) 






(4) 



with r*, r^^^ 



r^*^, and Mj denoting total width, partial 
width to hadrons, partial width to e+e^ pair, and mass of 
the resonance j, respectively. The total amplitude, which 
is the coherent sum of the three resonances, once squared, 
contains the interferences of the type ^T*Tj. Here (j>j is 
a phase associated to resonance j. If the resonances are 
quite broad, the interference effect will distort the shape 
of the resonances, the width might appear broader or 
narrower, and the position of the peak might be shifted 
as well. 

The total cross section of the resonances 
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FIG. 1; Total hadronic cross section in nb obtained as 
a{e'^e~ — >■ hadrons) — R ■ 86.85/s (s in GeV^) from R mea- 
surements of BES [111 . Iia ] and the fit to data. The curves 
show the best fit (identical for all the four solutions), the con- 
tribution from each resonance, as well as the interference term 
(one of the four solutions). 



and ^(4415) become wider than the previous published 
results where the interferences between resonances are 
neglected p^ . 

The fit results are presented in Table U there are four 
solutions found in the fit. It should be noted that the 
four solutions have identical x^, masses, and total widths 
for the resonances, but different partial widths to lepton 
pairs. One may suspect that the existence of four solu- 
tions is due to low precision of the measurements. We 
will show in the next session that multiple solution is a 
real effect. The improvement of the precision of the mea- 
surements can not change the fact that there are four 
solutions in fitting these data. 



III. MULTIPLE SOLUTIONS 
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Since what we actually obtain is the squared module of 
the amplitude, only two relative phases are relevant. 

By minimizing the y^ defined in Eq. ([Ij , we obtain the 
results as displayed in Fig.[T] Just as expected, the inter- 
ference effect changes the shape of each resonance signifi- 
cantly: the ■0(4040) becomes narrower while the ^'(4160) 



The existence of four solutions in the fit described 
above can be tested through the toy simulated data. The 
special steps for this approach is: 

1. The curve of the best fit in previous session (with 
background contribution subtracted) is used as the 
probability density function (PDF) to compute 
a set of experimental points at different energies 
{Npt = 100 in this work) within the range 3.75 GeV 
< V^ < 4.85 GeV. 



TABLE I: Four groups of solutions for the data fitting. Tfie 
four solutions have identical resonance masses (M) and total 
widths (Ft), but significant different leptonic partial widths 
(Fee) and the relative phases {(j>). The fit yields x^ ~ 91) 
A = (15.05 ± 0.59) nb, and B = (-1.64 ± 0.67) nb/GeV for 
all the solutions. 

Parameter V>(4040) ^/'(4160) ^(4415) 



M (MeV) 
Ft (MeV) 



4034 ± 6 
87 ±11 



4193 ± 7 
79 ±14 



4412 ± 15 
118 ±32 



Fii' (keV) 
4>^^' (radian) 



0.66 ±0.22 0.42 ±0.16 0.45 ±0.13 
(fixed) 2.7 ± 0.8 2.0 ± 0.9 



fW (keV) 0.72 ±0.24 0.73 ±0.18 0.60 ± 0.25 

<?!>(2) (radian) (fixed) 3.1 ±0.7 1.4 ± 1.2 

Fie' (keV) 1.28 ±0.45 0.62 ± 0.30 0.59 ± 0.20 

0^^' (radian) (fixed) 3.7 ± 0.4 3.8 ± 0.8 



TJT 



a(4) 



(keV) 1.41 ±0.12 1.10 ±0.15 0.78 ±0.17 
(radian) (fixed) 4.1 ±0.1 3.2 ± 0.3 



TABLE If: Fit results for four groups of solutions with the 
toy simulated data points. The definitions of the parameters 
are the same as in Table H) The fit yields x^ = 1-0 x 10"^ for 
all the solutions with the expectation of x^ ~ 0. 



Parameter 



V'(4040) 



V'(4160) 



V'(4415) 



M (MeV) 
Ft (MeV) 



4033.5 ± 0.3 
87.23 ± 0.49 



4192.8 ±0.3 
79.00 ± 0.53 



4412.4 ±0.4 
118.11 ±0.56 



'm' 



/,(!) 



=^7 



(keV) 
(radian) 



0.664 ± 0.005 
(fixed) 



0.417 ±0.004 
2.701 ±0.012 



0.454 ± 
2.002 ± 



0.003 
0.012 



A(2) 



^^137 



(keV) 
(radian) 



0.723 ± 0.006 
(fixed) 



0.731 ± 0.005 
3.051 ±0.001 



0.596 ± 
1.432 ± 



0.003 
0.014 



(keV) 
(f)'^^^ (radian) 



Fl!^(keV) 
(f)'''^^ (radian) 



1.283 ± 0.005 
(fixed) 



0.620 ± 0.006 
3.732 ± 0.006 



0.590 ± 
3.789 ± 



0.003 
0.013 



1.397 ± 0.006 
(fixed) 



1.087 ± 0.008 
4.082 ± 0.005 



0.774 ± 
3.218 ± 



0.003 
0.009 



2. A relative error of 1% and an absolute error of 
0.01 nb are added in quadrature and the total error 
is assigned to each data point. The inclusion of a 
0.01 nb absolute error is to weaken the chisqure- 
weight of points with small cross sections. 

Fit the simulated data with the similar x^ in Eq- (H} 
but with only resonance cross section included (i.e., 
cr*'"=(s) — a^'^'^-{s) in Eq. ([5])), the four groups of solu- 
tions obtained are summarized in Table |II] and displayed 
in Fig. [5] Comparison of two tables obviously indicates 
that the central values of the parameters are consistent 
with each other. 

From Table |TT] and Fig. [51 we can see that the largest 
Fee is more than twice the smallest value in the four so- 
lutions. We also notice that the first solution is identical 
to the BES published one [131 , and the second solution is 
about the same as the one listed as the best estimation 
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FIG. 2: Four groups of solutions obtained from the fit to the 
generated data. For clearness, the interference curves have 
been moved downward by 10 nb, the black dashed line at 
— 10 nb corresponds to zero cross section. 



of the partial widths of these states by the PDG [10| . 



IV. DISCUSSIONS 

As the Fee of the vector resonances are closely related 
to the nature of these states, the choice among the dis- 
tinctive solutions affects the classification of the charmo- 
nium and charmonium-like states observed in this energy 
region. 

The calculation of the e~^e~ partial widths of the S- 
wave charmonium states is well summarized recently in 
Table 2 of Ref. [16]. We can see clearly that many of 
the theoretical calculations give large ipiSS) and ■0(45') 
decay widths compared to the PDG values (about the 
same as the second solution listed in Table U above) , but 
the agreement with the third or the fourth solution is 
much better. 

The F(4260) was proposed to be the ^(45) state 
and the 7/;(4415) be ip{5S) in Refs. [3, [9l|, we can see 
that in this assignment, the calculated partial widths of 
V'(4040) = tPiSS) and ^(4415) = V^SS*) Q agree well 
with the fourth solution listed in Table H] 

Of course the possible mixing between S- and D-wstve 
states will change significantly the theoretical predic- 
tions of the partial widths of these states (iTII , and the 
QCD correction, which is not well handled [i6[, may also 
change the theoretical predictions significantly. So far, 
we have no concrete criteria to choose any one of the 



solutions as the physics one. 

It should be noticed that if the Y states are considered 
together with the excited tp states in fitting the i?-values, 
there could be even more solutions, the situation may 
become more complicated. 

We also notice that the existence of the multiple solu- 
tions is due to the inclusion of a free phase between two 
resonances, if these phases can be determined by other 
means (either theoretically or experimentally), it will be 
very helpful to know which solution corresponds to the 
real physics. 



V. SUMMARY 

Based on the R scan data, the resonance parameters 
of excited '0-family resonances are fitted. We found that 



there are four sets of solutions with exactly the same fit 
quality extracted from the experimental data, but the 
leptonic partial widths among different sets of solutions 
differentiate from each other significantly. New informa- 
tion is needed to determine which solution corresponds 
to the real physics. 
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